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Abstract 

It is shown that charged hadron multiphcity distributions in lepton-nucleon 
scattering are fairly well described by the modified negative binomial distri- 
bution in the energy range from 3 — 4 to ~ 220 GeV. The energy behaviour 
of the parameter k is similar to the dependence observed for e~^e~ annihi- 
lation. 
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In recent papers [1,2,3,4,5] negatively charged particle multiplicity distributions 
in e+e^ annihilation have been well described by the modified negative binomial 
distribution (MNBD) with the parameter A practically energy independent and 
the parameter k approaching the value 7. The aim of this letter is to show 
that the MNBD fairly well describes multiplicity distributions in lepton-nucleon 
interactions with similar behaviour of the parameters. 

Let us remind that the MNBD can be defined by the probability generating 
function 



where P„ is the probability to produce n particles. For integer k, positive r 
and negative A lying in the interval — 1 < A < the MNBD is positive definite 
and can be computed as a convolution of the Newton binomial (the numerator in 
(1)) and the negative binomial (the denominator in (1)), these two distributions 
are also positive definite. The probabihties P„ for the MNBD can be computed 
also using iteration relations given in the paper [4] or the formulae given in the 
paper [6]. 

In tables 1 and 2 the results of the MNBD fits with fixed A = —0.5 to the 
published negatively charged particle multiplicity distributions for the reactions 
^± + p ^ ^± + X+ [7,8,9] and /i+ + ^ + X+ [10] are given for different 
intervals of the invariant mass W of the hadronic system X and the momentum 
transfer squared Q^. The parameter r has been replaced in fits by the the mean 
multiplicity < n > using the constraint 



The data have been selected using two simple criteria: a) the number of measured 
points Up in the multiplicity distribution exceeds 3 and h) W > 3 GeV. For 
reaction the published uncorrected multiplicity distributions have been used, since 
the effect of corrections is rather small as can be seen from the comparison of the 
uncorrected multiplicity distributions with results of the negative binomial fits to 
the "unfolded" distributions [10,11]. 

One can see from tables 1 and 2 that the quality of fits is quite qood, the 
points with big /NDF values are mainly at small W and/or Q'^ . Big x^/NDF 
at low W can be explained by the rise of the parameter k from ~ 1 to ~ 3 in the 
wide W intervals, at higher W the relative variation of the parameter k becomes 
smaller and therefore its influence is not seen. It is necessary to note that for 
points with big /NDF acceptable fits are obtained with free A (not shown). 
The reasonable fits with free A are obtained also in the energy range 1^ < 3 GeV, 




(1) 



r = A+ <n> /k. 



(2) 
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where Up is usually below 4, i.e. the number of degrees of freedom is too small 
for two-parameter fits. 

In the table 3 the results of the MNBD fits to the multiphcity distributions 
in the reactions z/ + J9 — > /i~ + [12,13], u + n —>■ IJ,~ + X+ [12] and V + p ^ 
yU"*" + [14] are given. For these reactions the quality of fits is also quite qood, 
big /NDF again arc observed mainly for the W intervals near 3 — 4 GcV. 

As seen from tables 1 and 2 the dependence of the parameter k for ixp 
reactions in different W intervals is weak or even absent. The energy dependence 
of the parameter k for the studied lepton-hadron processes is compared in fig. 1 
with the energy dependence of the parameter k for e^er annihilation, the k values 
for e+e^ data were taken from the papers [4,5]. For clarity only the part of the 
values given in tables 1-3 is shown in this figure. The errors of the parameter k 
for the studied reactions (not given) are of the order ~ 1 and exceed the errors 
for statistically more significant e^e~ data. The k values for reactions follow 
the dependence observed for e+e~ annihilation. The k values for vp scattering 
are slightly below the ones for e+e~ annihilation and approach them from below 
with the energy increase whereas the k values for Vp scattering are slightly above 
the e+e~ data. The k values for ji^D^ scattering are shghtly below the k values 
for scattering. 

Recently charged particle multiplicity distributions for the reaction + p ^ 
+ X have been published by the HI Collaboration [15]. The results of the 
MNBD fits with the fixed k = 7, as suggested by the e'^e^ data, are given 
in the table 4 for the maximum available pseudorapidity interval 1 < 77* < 5. 
The /NDF values for MNBD fits are good, but they should be considered 
just indicative, since the full covariance matrix is not published and therefore 
the proper treatment of the correlations between measurements of the neighbour 
multiplicities is not possible. This contrasts presented earlier fits of the muon 
and ( ant i) neutrino data obtained mainly using bubble and streamer chambers, 
where these correlations are small. Nevertheless one can compare these values 
with the /NDF^ obtained for fits using another generalization of the negative 
binomial, so called Generalized Negative Binomial Distribution (GNBD) [16,17] 
(see table 4). The /NDF values for MNBD fits are significantly smaller than 
for GNBD fits. It is necessary also to note that MNBD fits are one-parameter 
fits, since only A is free parameter and r has been fixed using the measured mean 
multiplicity < rich >, whereas the GNBD fits are two-parameter fits. 

One can also estimate the parameter k from the ratio of the dispersion D — 
{< v? > — < n >^)^/^ to the mean multiplicity < n> using formula 



L>2 



1 



1 + 2A 

< n > ' 



(3) 



< n >2 
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for A = —0.5 it gives 



This formula has been apphed to the data given in the papers [18,19,20], 
where multiphcity distributions are not tabulated. The ratio of the dispersion 
Dch to the mean multiplicity for charged particle multiplicity distribution has 
been calculated from the results of the negative binomial fits (parameters k]sfB 
and < n >nb) and has been transformed to the ratio for negative charged particle 
multiplicity distribution. The final formula for the parameter k (called k to be 
distinguished from k for MNBD fits) is 

/ k^B < n >NB N/< n >NB ~\(l\^2 /r\ 

i^={l )( )) (5) 

kNB+ < n >NB < n >NB 

where q is the total charge of the hadronic system X in the studied reactions. 

The energy dependence of the parameter k shown in the fig. 2 is similar (if 
not coincides) with the energy dependence of the parameter k in the fig. 1. The k 
values for the HI data, shown in the figure were calculated from the values of the 
second normalized moment C2 =< n'^ > / < n given in [|15| for the interval 



1 < Tj* < 5 using the relation (< n > /-D)^ = 1/(C*2 — 1). The C2 for the charged 
particle multiplicity distribution are used, since at the HERA energies the effect 
of the total charge q of the system X (second factor in the formula 5) can be 
neglected. 

What is the physical meaning of the MNBD parameters k and A? One can 
assume that k is the number of sources of particle production at some stage of 
interaction, each source develops according to some branching process produc- 
ing intermediate neutral clusters (as proposed in [Q, |, |])0. The constancy of 
the parameter A indicates that this branching process is pure birth branching 
process and in this case —A represents the cluster decay probability into charged 
hadron pair (1 + A is the probabilty of cluster decay into pair of neutral hadrons). 



Following the statistical model proposed by Goulianos it was assumed in 
[|^, 1^, H that A should be equal -0.5. But this assumption is incorrect, since it 
predicts wrong ratio of the average neutral and charged pion multiplicities. In the 
Goulianos model this deficiency was corrected by "asymmetric" assumption that 
the neutral cluster decays into charged pion pair or transforms into one neutral 
hadron. The observed neutral/charged ratio is restored in the model with 27r° 
cluster decay mode if A = —2/3, this value is close to the value A = —0.76 used 
in the MNBD fits of the e+e" data. One should note that the value A = —2/3 is 
expected if neutral clusters have zero isotopic spin. The fits of the lepton-nucleon 



•^It is of interest to note that the parameter k from MNBD fits is rather close to the parameter 
c of the Levy parametrization for the KNO scahng function V'(-z) f^i [ill, ^ • The possible relation 
between Levy and MNBD parametrizations has been mentioned in 
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data with A equal to -2/3 or -0.76 give also acceptable fits but in general worse 
compared to the A = —0.5. 

In conclusion, it is shown that the MNBD with fixed parameter A = —0.5 
quite well describes lepton-nucleon multiplicity data in the wide energy range 
from 3 — 4 to ~ 220 GeV. The parameter k for lepton-nucleon interactions has 
the energy dependence similar to the energy dependence for e~^e' annihilation. 
The evidence, given in this paper, for the same asymptotic value k = 7 for e~^e~ 
annihilation and lepton-nucleon scattering has no direct theoretical explanation 
now. Future experimental and theoretical studies should show whether this co- 
incidence is accidental or reflects some common dynamical mechanism. 
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Figure Captions 



Fig.l The energy dependence of the parameter k obtained from the MNBD fits 
with fixed A = —0.5 to the negative charged particle multiphcity distributions in 
the lepton-nucleon interactions compared to the k values for e~^e~ annihilation. 
For clarity some points are slightly shifted from integer values. The histograms 
are drawn to guide the eye. 



[<n> /Df = l/{C2-1) 

compared to the k values for e^e~ annihilation. 
For the HI data the measurements in the pseudorapidity interval 1 < rj* < 5 are 
used. 



Fig. 2 The energy dependence of the parameter n = 
for the reactions e~^+p e++X [0, fi^+p fi~^+X 
and u + p 
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Table 1: Results of the modified negative binomial fits with fixed pa- 
rameter A = —0.5 to the negative charged particle multiplicity distributions 
in lepton-hadron reactions fi"^ + p ^ fi'^ + X [9], + p ^ fi~ + X [7] and 
fi+ + D2^ fi+ + X [10]. 



Experiment 


W (Gel/) 


Q\GeVf 


< n_ >(fit) 


k 


XVNDF 










0.24 


1.146±0.018 


2 


10.9/3 






2.8- 


=-3.8 


0.54 


1.017±0.018 


2 


7.0/3 










0.95 


1.071±0.024 


2 


7.0/3 


J.Ballam 






1.80 


0.782±0.076 


1 


13.1/3 


^M-p 






0.0675 


1.222±0.029 


3 


30.4/3 


16 GeV/c 


3.8- 


=-5.0 


0.25 


1.181±0.075 


2 


5.8/3 










0.56 


1.354±0.131 


2 


0.8/3 










1.38 


1.151±0.056 


3 


0.2/3 






4- 


=-6 




1.754±0.071 


4 


1.3/5 






6- 


r8 




2.067±0.086 


3 


3.3/5 


NA9 


i 


8^ 


-10 




2.339±0.037 


4 


9.4/6 




10- 


rl2 


4^200 


2.598±0.080 


5 


19.2/6 


280 GeV/c 


12- 


^14 




2.848±0.083 


5 


6.8/6 






14- 


=-16 




3.011±0.072 


6 


7.2/7 






16- 


=-18 




3.108±0.081 


6 


19.1/8 






18- 


r20 




3.206±0.068 


7 


39.6/8 






8^ 


-10 




4.984±0.108 


4 


4.9/14 






10^ 


-12 




5.474±0.095 


4 


17.6/15 






12-^ 


-14 




5.943 ±0.106 


4 


33.3/17 


E665 


1C| 


14- 


hl6 




6.525 ±0.115 


4 


21.4/19 




2 


16- 


=-18 


> 1.0 


6.787 ±0.112 


4 


36.4/21 


490 GeV/c 


18- 


=-20 




7.336 ±0.131 


4 


30.4/19 






20- 


=-22 




7.241 ±0.131 


5 


22.2/21 






22- 


=-25 




7.813 ±0.135 


4 


36.4/23 






25- 


=-30 




7.861 ±0.156 


4 


28.3/22 
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Table 2: The same as in Table 1 for reaction //+ + p — > + X at 14 GeV/c [8]. 



W (OpV) 




\ / 1/ — ^ yiLii J 


]r 




O.i 




1 0704-0 n7K 
i.U / ZztU.U ( 


9 


o.U/ 


Q K 
O.O 




i.ozo±U.ic)0 


Q 
O 


o.U/o 


Q 7 


U.o— U.o 


1 ooc^-un 117 

i .ZZOltU. i i ( 


O 

A 


O 7 /Q 
Z,. 1 1 6 


o.y 




i.ozOitU.ioo 


Z 


T A O /"i 

l^.Zj 6 


/I Q 




1 Q7o_i_n 1 1 (S 
i.o /oitU.iio 


4 


O.O/ 6 


o. i 




1 nno-un iTifi 

i.UUZitU.UOO 


o 
Z 


iO.o/ 


Q Q 
O.O 




1 no/i -Ln n7n 


o 
z 


U.i/o 


O.O 


U.Do 


1 noo_i_n n/^n 
i.Uzo±U.UDy 


o 
z 


7 O /Q 

/ .z/o 


C). I 




1 on7^n nc^c; 
i.zU < ±U.Udd 


Q 
O 


u.y/o 


o.y 




1 onQ^n n7n 

i.zUo±u.u /y 


o 
z 


z.o/o 


^. o 




1 Q«o 1 n 1 oi 
± .ouonzu . ±o 1 


9 


O.U/O 


3.1 




1.405±0.052 


2 


12.4/3 


3.3 




1.119±0.024 


3 


1.3/3 


3.5 


1.75 


1.061±0.038 


2 


10.0/3 


3.7 




1.150±0.075 


2 


6.3/3 


3.9 




1.141±0.180 


1 


1.5/3 


4.3 




1.328±0.121 


2 


3.7/3 




0.3^0.5 


1.042 ±0.037 


2 


25.3/3 




0.5^1.0 


1.002 ±0.033 


2 


17.4/3 


2.8^3.6 


1.0^1.5 


0.987 ±0.033 


2 


7.7/3 




1.5H-2.0 


1.027 ±0.037 


2 


11.3/3 




2.0^4.5 


1.054 ±0.038 


2 


1.4/3 




0.3^0.5 


1.296 ±0.074 


3 


23.7/3 




0.5^1.0 


1.358 ±0.074 


2 


10.2/3 


3.6H-4.7 


1.0^1.5 


1.163 ±0.052 


2 


4.5/3 




1.5^2.0 


1.279 ±0.105 


3 


1.4/3 




2.0^4.5 


1.261 ±0.147 


2 


4.8/3 
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Table 3: The same as in Table 1 for up ^ fi + X~^~^ [12,13], un ^ + [12] 
and up ^ fi~^ + X° [14] reactions. 



Experiment 


VV I UrC K I 


\ It — ^ I ilL I 










2.83—3.46 


0.766±0.022 


1 


32.6/3 






C).4D— 4.U 


1 1 QD-i-n no/i 
i.iC)U±U.U/4 


o 
Z 


Q n /c; 
o.U/0 






/I n • /I o 


i.2o ±U.030 


o 
Z 


1 o / 
2i.o/ 


up 




4. o— 0.66 


i.606±0.036 


Q 
O 


0. 1 /o 


D. Zieminska |[T^ 


O.DD— D. lb 


1 Q^n_i_n n /1 7 
i.ooU±U.U4^ 


Q 

6 


1 O /I /i^ 

iZ.4/0 






0. (0— ( .y4 


Z.U4oitU.Uuo 


Q 

o 


io.z/ 






7 0/1 • o /in 




/I 

4 


n 7 
y. / /D 






n /I n • 1 1 IS 

y.4y— ii.io 


Z.D4z±U.i/4 


A 

4 


1 Q 1 /fi 

io.i/0 






11 1 fi • 1 c; n 
ii.io— iO.U 


o 71 ri-un 1 c^i 
z. ( iUitU.ioi 


A 

4 


o.o/o 






O OO • O /I ^ 

2.00— 3.4d 


1 no^_i_n mo 
i.036±0.0i2 


o 
2 


O/I r\ 1 o 
24.U/3 






o /I . /I n 

3.46— 4. U 


1.185±0.020 


2 


O C O / /I 

25.2/4 






4.U— 4.0 


1 c;«n_|_n non 


Q 
O 


io.u/ o 


un 




4.8^5.66 


1.721±0.030 


3 


8.1/5 


D. Zieminska 


12 


1 


5.66^6.76 


1.990±0.048 


3 


3.7/5 






6.76^7.94 


2.217±0.054 


4 


1.7/5 






7.94^9.49 


2.476±0.076 


4 


3.4/6 






9.49^11.18 


2.004±0.036 


5 


9.8/6 






11.18^15.0 


3.128±0.187 


3 


5.8/6 






3.16-^4.47 


1.065 ±0.065 


2 


2.8/3 


up 




4.47^5.48 


1.279 ±0.085 


2 


2.8/4 


J. Chapman |]TI 


1 


5.48^7.07 


1.776 ±0.102 


3 


15.3/4 






7.07^10 


1.564 ±0.130 


2 


9.0/6 






3.0^3.5 


1.196 ±0.042 


2 


7.5/3 






3.5^4.0 


1.545 ±0.033 


3 


9.4/3 


up 




4.0^4.5 


1.864 ±0.069 


4 


3.7/3 


M. Derrick [14 


1 


4.5^5.0 


1.899 ±0.088 


4 


3.4/3 






5.0^7.5 


2.510 ±0.073 


5 


13.3/4 






7.5^10.0 


2.840 ±0.175 


8 


0.3/4 
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Table 4: The results of the MNBD fits with k = 1 to the HI data [15] with 
< rich > fixed at the published values for the pseudorapidity interval 1 < 77* < 5 
and for the interval 10 < < iqoO (Gev/c)^. In last column the x^/NDF 
for fits using Generalized Negative Binomial Distribution from the paper [17] are 
shown. 



W (Gev) 


-A 


XVNDF 


xVNDF(from [17]) 


80^115 


0.216 ±0.034 


2.7/18 


12/17 


115^150 


0.233 ±0.036 


10.6/21 


26/20 


150^185 


0.187 ±0.049 


4.6/22 


19/21 


185H-220 


0.264 ±0.057 


11.3/23 


16/22 



11 



10 



7 - 



6 - 



5 - 



4 - 



3 - 



2 - ^-B 



■A- 



-A---A- 



■Q- 



— a- 



— a- 



4 



r--A^T--A- 



O /:x^p ^ /i,^ + X, NA9, Arneodo et al. 

® /j.^ D 2 ^ /J.^ + X , E665, Adams et a 

□ i/^p ^ /x^ + X, Derrick et al. 

A i/^p ^^t" + X, Zieminska et ol. 



3 



4 5 6 7 8 9 10 



20 



GeV 



30 



12 



10 



A 
V 



7 - 



6 - 



5 - 



4 - 



3 - 



2 - 




-A- 

-A- 



-A- 



X 




O e"p^e^+X, HI , Aid et al. 

O fi'^p — >/x^ + X, NA9, Arneodo et al. 

□ i/^p ^/x^ + X, WA21 , Jones et al. 

■ y^p ^z^^ + X, WA25, Jongejons et al, 

A L/^p ^^a^ + X, WA21 , Jones et al. 

A i/^p — >/i." + X, WA25, Jongejans et a 
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